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ABSTRACT 


The  propane-water  system  was  studied  in  three  phase 
(hydrate  -  water  rich  liquid  -  gas)  equilibrium  in  the  region 
of  32 °F  to  42.3°F  and  25  psia  to  80  psia.  This  equilibrium 
was  terminated  at  the  four  phase  point  due  to  the  appearance 
of  the  propane  rich  liquid  phase  at  T  -  42.3°F  and  P  =  80  psia. 

The  carbon  dioxide  -  water  system  was  investigated  in 
the  same  three  phase  equilibrium  in  the  region  of  32 °F  to 
50.3°F  and  172  psia  to  648  psia.  The  quadruple  point  was 
observed  at  T  =  50.3°F  and  P  =  648  psia. 

The  ternary  system  co2~c3hq~h2°  was  studied  in  the 
three  phase  equilibrium  region  at  various  compositions.  The 
four  phase  locus  was  determined  and  a  maximum  temperature  point 
was  observed  in  this  locus  at  56.7  °F. 

The  potential  parameters  for  CO^  and  C-^Hg  molecules 
were  computed  by  fitting  the  experimental  dissociation  pressure 
at  0°C.  They  were  6  =  3.96  A  and  £/k  -  148  °I<  for  the  carbon 
dioxide  molecule,  while  d  =  4.80  A  and  £/k  -  278  °K  for  the 
propane  molecule.  Using  these  parameters,  Langmuir  constants 
were  calculated  as  a  function  of  temperature.  These  constants 
were  used  to  predict  initial  hydrate  formation  conditions  for 
CC^-H^O  and  C^Hg-E^O  systems.  The  average  deviations  between 
the  predicted  and  experimental  results  were  6%  and  4%  for  the 
CO^-H^O  and  CgH^-EUjO  systems  respectively. 
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INTRODUCTION 


A  gas  hydrate  is  an  inclusion  compound  in  which  two  or  more 

components  are  associated  without  ordinary  chemical  union,  but 

through  a  complete  enclosure  of  one  set  of  molecules  in  a  suitable 

structure  formed  by  another.  Such  compounds  were  classified  as 

(  35) 

Clathrates  by  Powell  .  They  are  characterized  by  a  cage-like 
structure  which  is  formed  by  water  molecules.  The  hydrate 
forming  molecules  are  enclosed  in  these  cavities. 

The  process  of  hydrate  formation  can  be  represented  as, 


R 


(g) 


+  n  HO  v 

2  (1) 


n  H  O, 

2  (s) 


where  a  hydrate  forming  gas  R  combines  with  n  molecules  of  water 
to  form  a  hydrate  R  •  n  H2O.  As  the  forces  stabilizing  the  hydrate 
are  purely  physical  in  nature,  a  gas  hydrate  is  not  a  stoichio¬ 
metric  compound.  The  experimental  value  of  n  varies  from  5.70 
to  17.0  depending  upon  the  hydrate  forming  molecule  and  the 
structure  of  hydrates  formed. 

The  type  of  molecules  which  may  be  bound  in  the  cavities 
is  very  diverse,  as  it  includes  O2 /  N2/  CH^,  C2H2,  CH3CF2CI  etc. 
and  even  the  noble  gases  A,  Kr,  Xe .  Once  such  a  molecular  complex 
is  formed,  it  may  persist  under  conditions  where  it  is  no  longer 
thermodynamically  stable.  An  example  is  the  complex  with  Argon 
which  can  be  kept  indefinitely  in  an  ordinary  bottle,  although  the 
equilibrium  pressure  of  argon  over  the  crystal  amounts  to  several 
atmospheres  at  room  temperature. 


2  - 


Faraday,  in  1810,  discovered  the  first  gas  hydrate  of 

chlorine.  Villard^^  and  de  Forcrand^^"'  ^ studied  the 

hydrate  forming  conditions  and  hydrate  formulas  of  hydrocarbor 

( 15) 

and  non-hydrocarbon  compounds.  Hammer schmidt ' s  publications 

confirmed  that  gas  hydrates  were  responsible  for  the  plugging 

of  natural  gas  transmission  lines.  His  work  provided  an  impetus 

to  research  in  the  area  of  process  development  to  inhibit  hydrate 

formation.  Methods  and  means  were  proposed  for  the  partial 

dehumidification  of  the  natural  gas.  This  technique  was 

successfully  used  to  prevent  condensation  of  water  vapors  in 

the  transmission  lines,  at  operating  pressures  and  temperatures 

which  were  also  favourable  for  hydrate  formation.  Various 

(4,  16) 

dehydrating  agents  were  proposed,  one  of  which,  diethylene 

glycol,  is  now  being  used  extensively. 

To  be  able  to  predict  the  initial  hydrate  formation  conditions 

of  a  gas  mixture  of  known  composition,  the  concept  of  vapor- 

solid  equilibrium  was  proposed  by  Wilcox  et  al^O)  a  This 

(2  5) 

method  proved  useful  for  engineering  calculations  ,  although 
it  could  not  account  for  the  effect  of  composition  on  the  vapor- 
solid  equilibrium  constants  of  multi-component  systems. 

Claussen^^  and  von  Stackelberg^^'  studied  the 

crystal  structures  of  gas  hydrates  through  the  use  of  the  x-ray 
diffraction  technique.  All  hydrates,  according  to  these  studies, 
crystallize  in  either  of  two  structures,  Structure  I  or  Structure 
II.  Dimensions  and  other  details  of  these  structures  are  given 


«• 
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in  Appendix  B.  Although  it  is  difficult  with  this  information 
to  predict  exactly  which  molecules  will  form  gas  hydrates,  the 
dimensions  of  the  cavities  may  restrict  the  number  of  gases 
which  can  form  hydrates.  Any  gas  whose  molecular  size  is 

o 

greater  than  6.7  A  cannot  form  hydrates.  Theoretically  there 
may  not  be  any  restriction  on  the  smaller  molecules.  However, 
it  is  found  that  smaller  molecules  like  hydrogen  and  helium  are 
not  retained  within  the  frame.  Gases  like  ammonia  and  halogens 
which  show  specific  chemical  reaction  with  water  molecules  do 
not  form  gas  hydrates . 

Dimensional  limitations  and  geometrical  shape  of  the 

cavities  in  Clathrate  compounds  suggested  uses  in  purification 

(3,  36)  (44) 

and  separation  processes  .  Schaeffer  has  reviewed 

extensively  the  applications  of  clathration  as  a  unit  process  in 


the  petroleum  refining  and  petro-chemical  industries.  Parent 


(34) 


has  evaluated  the  proposal  for  storage  of  natural  gas  as  hydrates. 

Very  recently,  gas  hydrates  are  gaining  importance  as  potential 

(2)  .  (27) 

means  for  desalination  of  sea-water  .  Propane  was  studied 
on  a  pilot-plant  scale  and  preliminary  process  design  calculations ^ 45^ 
were  made  to  compare  the  feasibility  of  this  process  with  other 
processes  available  for  demineralization  of  sea-water. 

In  order  to  study  the  effect  of  one  gas  on  the  hydrate 
formation  conditions  of  the  other  gas,  ternary  systems  A-B-H^O 
were  investigated  by  Deaton  and  Frost^®'  ^ ,  Reamer  et  al^4^, 
Otto^32^  and  Jhaveri^20^.  As  carbon  dioxide  and  propane  are 


4 


constituents  of  the  natural  gas,  it  was  decided  to  study  the 

initial  hydrate  formation  conditions  of  the  systems  C02-H20, 

C3H8-H2O  and  CC>2-C3Hg-H20 .  The  three  phase  HL^G  equilibrium 

of  the  binary  system,  propane-water,  has  been  investigated  by 

( 42 ) 

Miller  and  Strong  and  of  the  binary  system  C02-H20  has 

.  (  52 ) 

been  studied  by  Unruh  and  Katzv  1  . 

The  initial  hydrate  formation  conditions  of  certain 

ternary  systems  containing  either  C0„  or  C  H  have  been  reported 

3  o 

in  the  literature  for  CH^-CgHg-H^^ ,  CH4-C02-H20^  52 ^  and 

( 40) 

C3I*8”C3H6“**2<“>  *  The  temperatures  of  the  hydrate  formation 

are  lower  than  the  critical  temperature  of  CO^  and  in  the 

ternary  system  CC^-CgHg-H^O .  In  other  words,  the  hydrate 

forming  gases  will  appear  in  the  non-aqueous  liquid  phase  L2 

at  certain  pressures.,  Such  a  system  would,  therefore,  exhibit 

a  univariant  four  phase  (HL^L2G)  equilibrium. 

(54,  55,  56) 

Platteeuw  and  van  der  Waals  and  Barrer  and 

(  3) 

Stuart v  1  proposed  a  statistical  thermodynamic  approach  to 
predict  the  phase  behaviour  of  gas  hydrates  in  three  phase 
equilibria,  HIG  and  HL^G .  This  solid  solution  theory  has  been 
used  for  pressure-composition  calculations  of  the  ternary  systems 
C3H8-H2S-H20  and  C3H8-CH4-H20  in  the  three  phase  HIG  equilibrium. 
It  was  decided  in  this  work  to  extend  this  theory  to  the  study 
of  initial  hydrate  forming  conditions  of  the  system  C02~C3H8-H2° 
in  three  phase  equilibrium  HL^G. 


THEORETICAL  PROGRAM 


A.  Phase  Behaviour  of  Gas  Hydrates 

The  total  number  of  phases  present  at  the  equilibrium 
pressure  and  temperature  in  any  system  is  given  by  Gibb's 
Phase-rule  which  states  that, 

F  +  P  =  C  +  2 

where , 

F  =  number  of  degrees  of  freedom 
P  =  number  of  phases  present 
C  =  number  of  components 

From  this  relationship,  it  may  be  deduced^2  ^  that  the 

addition  of  a  component  to  a  system  adds  a  degree  of  freedom  to 

all  features  represented  by  the  system.  In  other  words,  areas 

will  generate  volumes,  lines  will  generate  surfaces  and  points 

will  generate  lines.  These  features  will  have  one  extra 

degree  of  freedom,  but  they  will  not  be  changed  in  designation. 

That  is,  a  point  representing  vapor-solid-liquid  equilibrium 

for  a  one  component  system  will  generate  a  line  in  space  on  the 

addition  of  a  second  component,  but  this  line  will  represent  the 

equilibrium  between  the  same  phases,  namely  vapor-solid-liquid. 

Termination  of  a  particular  type  of  phase  equilibrium  may 

be  brought  about  by  the  addition  of  a  phase,  by  the  loss  of  a 

component,  or  by  attainment  of  critical  conditions.  The 

appearance  of  a  new  phase  decreases  the  degree  of  freedom  by  one. 

(  32 )  (  34) 

In  light  of  these  general  rules,  Otto  and  Parent  have 
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described  the  essential  features  of  the  phase  diagrams  of  binary 
and  ternary  systems . 

The  phase  behaviour  of  ternary  systems  such  as  A-B-H^O 

can  be  classified  into  two  types: 

(i)  Systems  in  which  hydrates  of  A-B--H2O  grow  in  one  structure. 

Systems  of  this  type  are  illustrated  in  Figure  1. 

They  exhibit  complete  miscibility  of  the  solid  solutions, 

that  is,  hydrates  of  A-H20  and  B~H20  are  formed  in  either  of 

the  structures.  The  ternary  system  CH^-r^S-BLjO  is  an  example 

of  this  type  of  phase  behaviour.  This  particular  ternary 

(  31) 

mixture,  as  studied  by  Noaker  and  Katz  ,  showed  a  continuous 
range  of  mixed  hydrates  of  Structure  I.  In  this  system,  since 
the  temperatures  of  hydrate  formation  were  higher  than  the 
critical  temperature  of  CH4,  a  methane-rich  liquid  phase  did 
not  appear  at  any  pressure.  On  the  other  hand,  in  a  system  like 
CO2-H2 S-H20  both  of  the  hydrate  forming  components  will  appear 
in  the  non-aqueous  liquid  phase,  L2 .  As  shown  in  Figure  1, 
hydrates  in  the  binary  systems  C02-H20  and  H2S-H20  grow  in 
Structure  I.  Lines  joining  two  pairs  of  quadruple  points, 
one  at  the  freezing  point  (HjIL]_G)  and  the  other  above  the 
freezing  point  of  water  (HjGL^L2),  describe  the  univariant 
four  phase  equilibria.  Due  to  lack  of  experimental  data,  the 
exact  shape  of  these  four-phase  locii  is  not  known. 

An  example  of  complete  miscibility  in  Structure  II 


hydrates  is  the  ternary  system  propane-propylene-water. 

(40) 

system  has  been  studied  by  Reamer 
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FIG-  I  P-T  PROJECTION  OF  PHASE  BEHAVIOUR  OF  THE  SYSTEM 
H2S-C02-H20  IN  WHICH  THE  BINARY  SYSTEM  H2S~H20 
AND  CO2-H2O  FORM  STRUCTURE- 1  HYDRATES- 
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(ii)  Systems  in  which  hydrates  of  A-B-H20  grow  in  both 
structures. 

These  are  illustrated  in  Figure  2,  and  are  generally 
known  as  double  hydrates,  characterized  by  the  partial 
miscibility  of  solid  solutions.  In  other  words,  hydrates  of 
the  binary  system  A-H2O  may  grow  in  Structure  I  and  those  of  the 
system  B-H20  may  grow  in  Structure  II  or  vice  versa.  The 
ternary  system  C02-C3Hq~H20  is  an  example  of  this  type  of  phase 
behaviour.  The  binary  systems  CO2-H2O  and  C3Hg-H20  form 
hydrates  of  Structure  1  and  Structure  II  respectively. 

When  propane  is  added  to  carbon  dioxide  at  the  quadruple 
point  (HjL^IG)  ,  the  pressure  will  increase  until  the  mixed 
hydrate  of  Structure  II  is  formed.  Thus  the  four  phase 
equilibrium  (HtIL-iG)  for  the  three  component  system  is 
terminated  at  the  quintuple  point  (HjHj jL-, L2G)  due  to  the 

** f 

appearance  of  a  new  hydrate  phase  The  three-phase 

equilibrium  branches,  (HjIG)^  and  (HjLiG)^  therefore  lie 
below  the  four  phase  branches  (H^Hj^IG)^  ^  and  (HjHjjL^G) A  B 
respectively . 

Owing  to  the  presence  of  a  solid  phase  Hj,  the  branches 

(HiL2G)a  and  (L-jI^G)  are  shown  below  the  vapor  pressure 

curve  (LoG)  for  pure  carbon  dioxide.  Since  water  and  propane 
A 

are  less  volatile  than  carbon  dioxide,  the  four  phase 
equilibrium  branch  (HjHjjL^G)  lies  below  (H^L^G)^.  As  a 
result,  the  quintuple  point  (H^H^L-^I^G)  A  B  i s  located  at  a 
lower  pressure  and  lower  temperature  than  the  quadruple  point 
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FIG-2  PHASE  BEHAVIOUR  OF  THE  TERNARY  SYSTEM  C02-C3H8-H20 
IN  WHICH  THE  SYSTEM  C02-H20  FORM  STRUCTURE  I  AND 
THE  SYSTEM  C3H8- H20  FORM  STRUCTURE  U  HYDRATES 
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As  it  is  not  possible  to  detect  the  presence  of 

Hj.  and  H  phases  individually,  there  is  no  experimental 

evidence  available  to  support  these  deductions. 

It  is  interesting  to  note  that  the  four  phase  equilibrium 

(HjjL^I^G)^  B  starts  from  the  .quadruple  point  (HjjL^I^G)^  of 

propane.  This  locus  can  be  terminated  only  at  the  quintuple 

point  (H-j-H  L  L  G)  where  a  new  phase  H  appears.  A  maximum 
1  II  1  2  A,B  I 

temperature  point  is  observed  on  the  four  phase  (H  LnL  G) 

XX  X  z  A  ^  x3 

equilibrium  curve.  Beyond  this  temperature,  hydrates  do  not 
form  in  Structure  II. 

Co-existence  of  the  H_  and  H  „  phases  was  confirmed  by 

I  II 

(56) 

Platteeuw  and  van  der  Waals  when  they  studied  pressure- 

composition  diagrams  of  the  CH^-CgHg-H^O  and  H^S-C^Hg-H^O 

systems.  An  interesting  exception  to  these  discussions  is 

provided  by  the  system  H  S-CH  CHF9-H  0  in  which  the  two 

z  3  z  2 

binary  systems  H  S-H  0  and  CH  CHF  -H  0  form  hydrates  of 

2  2  3  £  l 

Structure  I  only,  but  the  ternary  system  forms  hydrates  of 
Structure  II. 


B.  Solid  Solution  Theory  of  Gas  Hydrates 

A  common  feature  of  all  gas  hydrates  is  a  host  lattice  of 
water  molecules,  which  by  itself  is  thermodynamically  unstable, 
but  is  stabilized  by  inclusion  of  the  second  component.  The 
forces  binding  this  component  may  be  similar  to  the  intermolecular 
forces  in  liquids ( 54 \  It  seems  natural  therefore,  to  regard 
them  as  a  solid  solution  of  the  second  component  in  the  host 


lattice . 
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Platteeuw  and  van  der  Waals^3^  derived,  through  the  use 
of  statistical  thermodynamics,  a  set  of  equations  which 
describe  the  phase  behaviour  of  Clathrates.  The  basic  equations 
are, 


=  M.P  +  RT  V  In  (1  -  £  yk.) 

k 

i 


1  +  ^  CJi  Pj 

J 


C 

ki 


1 


k  T 


hki  (T, V) 
0k  (T) 


hki  ^v)-xk 


(1) 


(2) 


(3) 


Equation  (1)  is  a  generalized  form  of  Raoult's  law  when 

solute-solvent  intereaction  is  neglected.  The  total  chemical 

potential  |j.q  is  made  of  contributions  due  to  the  empty  lattice, 

p 

represented  by  p,  and  that  due  to  encaged  molecules,  represented 
by  the  summation  term  in  equation  (1). 

Equation  (2)  is  equivalent  to  Langmuir's  isotherm  for 
localized  adsorption  without  interaction  between  the  adsorbed 
molecules.  The  probability  of  finding  a  molecule  k  in  cavity 
of  i  type  is  given  by  yk^.  The  equilibrium  pressure  or 
dissociation  pressure  p  can  be  replaced  by  a  corresponding 

1C 

fugacity  term  if  the  system  behaves  non-ideally. 

Equation  ( 3)  defines  the  Langmuir  constant  which  is 

a  function  of  temperature  only  for  a  given  molecule.  The 


. 


' 


12 


calculations  of  Langmuir  constants  are  illustrated  in  Appendix 

C. 

In  order  to  calculate  the  partition  function  h,  .  (T,V) # 

(56) 

Platteeuw  and  van  der  Waals  applied  the  Lennard-Jones- 

Devonshire's  cell  theory  which  is  appropriate  in  this  case 
because  the  gas  molecule  is  moving  in  an  almost  spherical  cage 
made  up  of  water  molecules.  They  calculated  the  force  field 
by  using  L-J-D  12/6  potential  model.  The  agreement  between 
calculated  and  observed  values  of  the  dissociation  pressure  is 
quite  good  for  the  monoatomic  gases  like  A ,  Kr,  etc.  and  to  a 
lesser  extent  for  the  quasi-spherical  molecules  like  CH^,  SF^ . 
However,  the  agreement  becomes  very  unsatisfactory  in  the  case 
of  non-spherical  molecules  like  CC^ ,  C2H6  an<^  N2  * 

This  poor  agreement  may  be  due  to 

(i)  hindered  rotation  of  the  rod-like  molecules  in  their 
cavities , 

(ii)  failure  of  the  central  force  field  approximation  to 
describe  the  interaction  between  a  solute  molecule  and  a 
molecule  in  the  host  lattice, 

(iii)  the  chemical  potential  of  the  lattice  may  not  be 
independent  of  the  solute  occupation, 

(iv)  the  shape  and  size  of  encaged  molecules  cannot  be 
neglected  when  evaluating  overall  cavity  potential. 

( 29) 

The  third  and  fourth  reasons  were  investigated  by  McKoyv 
To  a  good  degree  of  approximation,  it  was  confirmed  by  McKoy, 
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that  the  chemical  potential  of  the  lattice  is  independent  of 

the  solute  occupation  of  the  cavity .  This  shows  that  the 

possible  lattice  distortions  appear  to  have  a  very  small  effect 

on  the  thermodynamic  properties  of  the  lattice.  Improved 

(29) 

dissociation  pressures  were  obtained  by  accounting  for  the 

effect  of  the  shape  and  size  of  the  molecule. 

To  obtain  the  coefficients  for  the  Virial  equation  of 
state,  various  potential  functions  are  used.  These  potential 
functions  are  based  on  a  particular  model  proposed  for  the  gas 
molecule.  Realistic  potential  functions  must  be  used  to  fit 
the  experimental  data  over  large  temperature  ranges.  In  addition, 
more  realistic  potential  functions  were  proposed  to  obtain 
force  constants,  (  £/k)  and  (  <5 )  .  These  constants  could  be 
used  to  predict  the  behaviour  of  the  gas. 

The  Lennard- Jones  (6-12)  potential  has  been  widely  used 
for  the  calculation  of  the  properties  of  the  matter  in  the 
gaseous,  liquid  and  solid  states.  This  potential  function 
is  defined  in  terms  of  the  force  constants.  These  constants 
may  be  determined  by  analysis  of  the  second  virial  coefficient 
data.  McKoy^29^  found  that  the  dissociation  pressures  of 
gas  hydrates  were  more  sensitive  to  the  parameters  of  the 
intermolecular  potentials  than  were  second  virial  coefficient 
data.  The  values  of  these  parameters,  (  £  A)  and  (  6),  found 
in  the  literature  were  very  different  even  for  the  rare  gases 
which  have  been  most  studied. 
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For  equilibrium  between  hydrate  and  the  stable  crystalline 
a- form  (ice)  of  lattice,  we  have, 


Let  A  (a  (T,P) 


(4) 

(5) 


Since  the  a  and  forms  represent  the  pure  crystalline 
states  of  the  lattice  former,  the  difference  in  chemical 
potential  (A  \x)  is  a  function  of  pressure  and  temperature 
alone . 

Substituting  equation  (1)  and  (4)  into  equation  (5), 
we  obtain. 


A[jl 

RT 


(6) 


Substituting  for  y^  .  in  equation  (6)  from  equation  (2),  we 
get. 


A|a  =  RT  £  \>i  i"  u  +  E  Cki  Pk)  (7) 

i 

Equation  (4)  is  an  essential  condition  for  three  phase 
equilibrium,  HIG,  below  the  freezing  point  of  water,,  At 
temperatures  above  the  freezing  point  however,  the  stable 
liquid  form  L1 (water)  exists  in  equilibrium  with  other  phases. 
For  equilibrium  between  water  and  hydrate,  we  have, 

M-q  =  4q 


(8) 


-  15  - 


(9) 


Since  some  gas  is  dissolved  in  the  liquid  phase  L^,  the 
chemical  potential  difference  in  equation  (9)  is  now  a  function 
of  pressure,  temperature  and  composition. 

Assuming  that  the  water  in  contact  with  the  hydrate 
behaves  like  an  ideal  solution,  we  have. 


=  ( P-Q )  +  RT  In  xQ  (10) 


(10) 


where  jj,q  is  a  function  of  pressure,  temperature  and  composition 
while  (p,g)  is  a  function  of  temperature  and  pressure  only. 


Substituting  equation  (10)  in  (9), 


Let  A|i'  =  M^q  “  (M-q) 


Then, 


(ID 


A (i  =  Ap, '  -  RT  In  Xq 

where  Ap.'  is  a  function  of  P  and  T  only. 
Substituting  equation  (7)  in  (11),  we  get, 


(12) 


+  In  x, 


Q 


i 


At  constant  temperature,  the  effect  of  pressure  on 
chemical  potential  difference  is  given  by, 

(dAp’/dP)T  = 


AV 


(13) 


■ 
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Integrating  equation  (13)  at  constant  T, 

A(X «  =  Am.'  +  AV'  (P  -  Pn)  (14) 

o  u 

Equation  (14)  enables  one  to  calculate  Ap,'  values  at 
any  pressure  P  and  for  any  component  with  respect  to  Aji^  and 
PQ  values  of  a  reference  component  at  the  same  temperature. 

It  is  assumed  that  AV '  is  independent  of  pressure.  The  values 
of  AV'  for  Structure  I  and  Structure  II  hydrates  are  reported 
in  the  literature  ^  ^ ^  .  The  quantities  Ap,^  and  P^  for  a 
reference  component  are  computed  by  using  equation  (12) . 
Methane  and  sulfur  hexafluoride  were  chosen , as  the  reference 
components  for  Structure  I  and  Structure  II  hydrates 
respectively . 

Using  this  approach,  Marshall,  Kobayashi  and  Shozaburo 
(26,  27,  46,  47) 

investigated  A-N2-H2O  and  CH^-A-E^O  systems 
in  three  phase  (HL^G)  equilibrium.  Their  theoretical 
calculations  together  with  the  experimental  results  indicated 
that  the  solid  solution  theory  could  be  satisfactorily  used 
for  the  hydrates  of  non-spherical  molecules  in  the  three  phase 
hydrate  -  water  rich  liquid  -  gas,  equilibrium. 
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EXPERIMENTAL  PROGRAM 


A.  Equipment 

The  equipment  used  for  this  work  was  essentially  the  same 
as  that  used  by  earlier  workers,  Snell and  Jhaveri^20^  for 
their  work  on  gas  hydrates.  The  major  components  of  this 
equipment  are  shown  in  the  schematic  diagram,  Figure  3. 

The  equilibrium  cell  was  made  from  a  Penberthy  liquid 
level  gauge  capable  of  withstanding  6000  psia  at  250°F  with  a 
test  pressure  of  9000  psia.  The  cell  body  was  made  of  316 
type  stainless  steel  and  had  glass  windows,  one  on  each  side. 
These  glass  windows,  manufactured  by  Corning  Glass  Works,  were 
rated  at  9000  psia  at  100°F. 

The  cell  was,  as  shown  in  Figure  3,  jacketted  by  a  lucite 
shell  of  3/8"  thickness  and  8"  O.D.  The  top  and  bottom  of  the 
cell  were  connected  to  the  outside  of  the  jacket  through 
rotary  seals,  one  on  either  side.  These  seals  rest  on  bearings 
and  prevent  leakage  of  the  cell  contents  when  the  cell  is 
rocking.  Gaskets  made  from  1/8"  thick  Teflon  sheet  were  used 
to  seal  the  glass-to-metal  connections.  The  high  pressure 
tubings,  valves  and  other  connections  were  supplied  by 
Autoclave  Co.  and  were  rated  at  10,000  psia  pressure.  The  cell 
could  be  rocked  through  an  arc  of  30-35  degrees  by  an  eccentric 
driven  by  a  1/6  H.P.  motor. 
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fig-3  schematic  diagram  of  the  equipment 
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B.  Pressure  Measurement  and  Generation 

A  hand  driven,  positive  displacement  mercury  pump  of 
265  cc  displacement  volume  and  rated  at  8000  psia  was  used  to 
change  the  pressure  in  the  system.  Chemically  purified  and  once 
distilled  mercury  was  used  as  the  confining  fluid.  The  volume 
displaced  could  be  read  from  the  graduated  scale  on  the  pump. 

Pressure  was  measured  by  a  0-1000  psi  Heise  gauge  having 
a  minimum  scale  reading  of  2  psi.  The  pressure  calibration 
of  this  gauge  was  checked  with  a  standard  dead  weight  tester. 

C .  Temperature  Measurement  and  Control 

To  facilitate  visual  observations,  a  transparent  fluid 
with  low  freezing  point  was  chosen  for  heat  exchange  purposes. 
This  fluid  was  a  petroleum  based  solvent  having  a  flash  point 
of  110 °F„  The  cell  contents  were  maintained  at  the  desired 
temperature  by  circulating  the  bath  fluid  around  the  cell 
through  the  lucite  shell .  The  temperature  of  the  bath  fluid 
was  controlled  by  a  bimetallic  strip  controller  supplied  by 
American  Instrument  Company. 

The  temperature  of  the  cell  contents  was  measured  by  an 
Iron-Con stan tan  thermocouple,  inserted  through  the  side  of  the 
stainless  steel  body.  A  potentiometer  was  used  to  measure  the 
voltage  difference.  The  calibration  curve  for  the  thermocouple 
is  shown  in  Appendix  A. 
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D .  Materials 

The  mercury  used  for  this  work  was  chemically  purified 
and  distilled  once.  The  water  used  for  the  experimental 
work  was  deaerated  and  doubly  distilled.  The  carbon  dioxide 
and  propane  used  for  hydrate  formation  experiments  were  99.8% 
and  99.5%  pure  respectively. 

E .  Experimental  Procedures 

The  experimental  procedures  included  preparation  of  the 
gas  mixtures,  observation  of  the  initial  hydrate  formation 
conditions,  and  analysis  of  the  gas  samples. 

(i)  Preparation  and  Analysis  of  Gas  Mixtures 

\ 

A  mixture  of  the  desired  composition  of  carbon 

dioxide  and  propane  was  prepared  by  adjusting  the  pressures 

of  the  gas  cylinders.  From  the  compressibility  data  for  the 

(1  21  39) 

individual  components  and  their  mixture  '  '  ,  these 

pressures  were  calculated  since  the  volume  of  the  sample  bottle 
and  the  temperature  were  known.  These  mixtures  were  then 
analyzed  chromatographic ally.  The  chromatographic  analysis 
is  described  in  Appendix  A. 

The  lines  and  sample  bottle  D  (See  Figure  3)  up  to  valve 
9  were  evacuated.  Propane  was  then  bled  through  the  system  and 
evacuated  again.  Propane  was  taken  first  into  the  sample 
bottle  D  until  the  gauge  on  the  propane  cylinder  indicated 
the  desired  pressure.  Carbon  dioxide  was  then  added  to  propane 


' '  vr 
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until  the  desired  pressure  of  the  mixture  was  read  on  the 
gauge.  Propane  was  charged  first  because  the  propane  cylinder 
was  at  a  lower  pressure  than  the  carbon  dioxide  cylinder. 

Valves  7  and  8  were  then  closed.  A  sample  of  this  mixture 
was  then  taken  from  the  sampling  point  and  analyzed  on  the  gas 
chromatograph . 

(ii)  Initial  Hydrate  Formation  Experiments 

Above  the  freezing  point  of  water,  for  three 
phase  equilibrium  HL^G ,  of  the  ternary  system  C02-C3Hg-H20 , 
two  degrees  of  freedom  are  available  to  the  system.  This 
equilibrium  will,  therefore,  generate  a  surface  on  the  P-T 
diagram.  Such  an  equilibrium  with  a  constant  gas  composition 
will  generate  a  line  on  this  surface.  Thus  for  a  fixed  composition 
parameter  of  a  ternary  system,  there  exists  only  one  combination 
of  pressure  and  temperature  when  hydrates  will  be  in 
equilibrium  with  the  gas  phase  and  aqueous  phase. 

For  a  given  gas  composition  one  can  calculate  hydrate 

formation  conditions  from  the  vapor-solid  equilibrium  constants 

(21) 

for  propane  and  carbon  dioxide  .  The  contents  of  the  cell 
are  then  kept  at  a  pressure  slightly  above  this  predetermined 
pressure.  The  system  is  then  supercooled  by  10-12°F.  The 
equilibrium  cell  is  then  rocked  till  the  hydrates  are  formed. 
Hydrate  formation  was  accompanied  by  a  drop  in  pressure.  The 
temperature  was  then  raised  very  slowly  until  a  few  crystals 
remained  on  the  glass  window.  Due  to  dissociation  of  the  hydrate 


. 
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crystals,  pressure  increased  with  an  increase  in  temperature* 

The  crystals  on  the  glass  window  were  examined  under  a  light 
source  through  a  magnifying  glass*  When  these  crystals  started 
decomposing,  the  temperature  was  lowered  by  1-2°F  to 
recrystallize  the  hydrates*  The  temperature  was  again  increased 
very  slowly.  The  temperature  and  pressure,  when  these  few  new 
crystals  dissociated,  were  noted*  At  this  time  a  sample  of 
the  equilibrium  gas  phase  was  taken  in  the  gas  burette  and 
analyzed. 

The  critical  temperatures  of  carbon  dioxide  and  propane 
are  87C8°F  and  206. 0°F  respectively.  Since  the  experimental 
temperatures  were  lower  than  these  critical  temperatures, 
these  components  will  liquify  in  the  phase  at  certain 
pressures*  To  obtain  these  conditions,  the  contents  of 
the  cell  were  compressed  at  about  75°F  until  L2  phase 
appeared.  Depending  upon  the  composition  of  the  feed  gas 
mixture,  the  composition  of  phase  will  vary  from  100% 
propane  to  0%  propane.  The  temperature  was  then  decreased 
and  the  equilibrium  cell  was  rocked  until  a  new  hydrate 
phase  appeared.  At  these  conditions  four  phases, 
were  present. 

A  gas  sample  was  prepared  as  described  in  part  (i) . 

Valve  2  was  then  opened  and  some  mercury  was  pumped  in  until  it 
appeared  in  the  bottom  of  the  cell.  The  valves  2,  4,  6,  7 
were  closed  and  the  system  evacuated  for  one  hour.  Valve  5  was 
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then  closed  and  about  20  ccs  of  distilled  water  were  taken  into 
the  cell  at  the  top  of  valve  1.  The  valves  6,  1  were  closed 
and  some  gas  mixture  was  bled  out  at  valve  4.  The  gas  mixture 
was  then  fed  into  the  cell.  The  valves  8  and  1  were  closed 
and  the  cell  connected  to  pressure  gauge  at  valve  2 .  The 
mercury  pump  E  was  then  manipulated  until  the  desired 
pressure  was  read  on  the  gauge.  The  temperature  was  set  on  the 
temperature  bath.  The  equilibrium  cell  was  then  rocked  until 
the  hydrates  were  formed.  The  crystals  on  the  surface  were 
spread  on  the  glass  surface  by  tilting  the  cell.  The  temperature 
was  increased  slowly  and  these  crystals  on  the  glass  window 
examined  carefully  until  they  just  started  to  decompose.  At 
this  point,  temperature  was  decreased  slowly  until  some  hydrates 
recrystallized  on  the  glass  surface.  These  crystals  were  then 
watched  carefully  through  the  magnifying  glass  when  the 
temperature  was  being  raised  very  slowly.  When  they  started 
to  decompose,  pressure  and  temperature  readings  were  noted. 

At  this  point,  a  sample  of  the  gas  phase  was  taken  into  a  gas 
burette.  The  lines  leading  to  the  sampling  point  were 
evacuated  prior  to  sampling  the  gas  mixture.  Care  was 
exercised  to  see  that  the  pressure  did  not  fluctuate  during  the 
sampling  procedure .  The  gas  sample  was  allowed  to  stand  in 
the  burette  for  about  30  minutes  and  then  analyzed  on  the 
gas  chromatograph. 


. 


M  i'iiV 


mew  eDfiv-. *4/.i  s  o.W  no  els;  ::  >  «-v  [T  , sittf  Bl&yrf  t;,:J 


24 


The  system  pressure  was  lowered  and  the  above  mentioned 
procedure  repeated „  For  each  mixture,  it  was  found  convenient 
to  start  with  the  quadruple  point  and  then  reduce  the  pressure 
for  additional  runs.  Care  was  also  taken  to  leave  some  crystals 
from  the  previous  run  to  serve  as  seeding  crystals  for  the 
next  run.  This  facilitated  hydrate  formation.  The  water  was 
changed  for  each  gas  mixture, 

F 0  Experimental  Difficulties 

Since  hydrate  formation  is  a  phenomenon  similar  to 
crystallization,  it  was  observed  that  until  crystallization  has 
commenced,  it  is  difficult  to  get  hydrates  to  form  even  at 
pressures  200  to  300  psi  higher  than  they  will  form  when 
crystals  are  present.  It  was  also  noticed  that  alhtough  most 
of  the  hydrates  decompose  readily,  a  small  amount  of  hydrates 
persist  for  a  considerable  time.  The  crystals  which  fail  to 
decompose,  are  either  submerged  in  the  water  or  completely 
covered  by  a  water  film.  Particularly  in  the  case  of  carbon 
dioxide  it  was  found  difficult  to  initiate  hydrate  formation. 

In  this  case,  therefore,  the  subcooled  contents  of  the  cell 
were  further  cooled  by  sudden  expansion  of  the  charge. 
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EXPERIMENTAL  RESULTS 


A.  Experimental  Accuracy  of  the  Results 

Three  variables  measured  in  this  work  were  pressure, 
temperature  and  composition 0  As  the  method  used  for  determining 
initial  hydrate  formation  conditions  was  visual,  the  results 
might  involve  some  personal  errors  also0  These  personal 
errors  were  minimized  by  repeating  the  procedure  in  the 
neighbourhood  of  hydrate  formation  conditions „  The  use  of  a 
magnifying  glass  also  facilitated  visual  observations 0 

The  errors  in  the  temperature  measurements  were  limited  to 
the  errors  of  the  sensing  device  and  the  measuring  device.  The 
iron-constantan  thermocouple  used  for  this  work  was 
calibrated  with  a  standard  mercury-in-glass  thermometer  at 
the  ice-point  and  room  temperature  of  72°F.  A  linear  correction 
was  applied  to  the  voltages  and  corresponding  temperatures 
read  from  the  potentiometer  conversion  tables.  The 
potentiometer  could  read  voltage  difference  equivalent  to 
±0 . 15°F • 

The  pressure  gauge  readings  were  corrected  for  level 
differences.  The  minimum  gauge  reading  was  2  psi. 

The  chromatographic  results  were  normalized  on  the 
prorata  basis.  The  observed  values  of  mole  percentages  of 
each  component  added  to  (100±3)%. 
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B.  Propane-Water  System 

Experimental  results  for  the  initial  hydrate  formation 
conditions  for  the  propane-water  system  are  tabulated  in 
Table  I  and  Table  II.  The  results  are  plotted  in  Figure  4. 
Investigations  of  other  workers  at  similar  conditions  are 
reported  in  Table  IV  and  compared  with  this  work  in  Figure  4. 

The  present  measurements  are  in  good  agreement  with  those  of 
other  investigators. 

The  range  of  pressure  (25-80  psia)  and  temperature 
( 32-42. 3°F)  in  which  HL^G  equilibrium  for  CgH^-H^O  system 
persists  is  quite  narrow.  The  three  phase  HL-^G  equilibrium 
is  terminated  at  the  quadruple  point,  observed  at  T  =  42.3°F 
and  P  =  80.0  psia.  The  vapor  pressure  curve  for  pure  propane 
is  shown  in  Figure  4,  to  give  an  idea  of  the  approximate 
position  of  the  HL2G  and  L]JL2G  phase  equilibria. 

C .  Carbon  dioxide  -  Water  system 

Villard^53^  was  the  first  investigator  to  work  on  this 
binary  system  in  three  phase  HL^G  equilibrium.  Later  Unruh 
and  Katz  ^ 52  ^  published  their  results  for  the  system  C02-CH4-H20 
which  included  their  experimental  observations  for  the  binary 
C02-H20  system. 

Experimental  results  of  this  work  are  tabulated  in 
Table  I  and  Table  II,  while  those  of  other  workers  are  reported 
in  Table  III.  These  are  plotted  on  a  semi-log  graph  in 
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Figure  5.  The  results  of  this  work  are  in  satisfactory 
agreement  with  Unruh 1 s  results,  but  Villard's  observations 

i 

reported  in  1894  do  not  agree  with  the  more  recent  data. 

The  three  phase  HL-^G  equilibrium  is  terminated  at  the 
quadruple  point  found  in  this  work  at  50.3°F  and  648  psia. 
The  vapor  pressure  curve  of  C02  is  shown  on  Figure  5  to 
indicate  the  location  of  HL2G  and  L^L2G  phase  equilibria 
of  the  binary  CO2-H2O  system. 

D.  Propane  -  Carbon  dioxide  -  Water  system 

The  initial  hydrate  formation  conditions  of  ternary 
systems  containing  either  C02  or  are  reported  by  the 

following  workers; 


Authors 

System 

Reference 

1. 

Deaton  and 

[  Frost 

ch4-c3h8-h2o 

(10) 

2. 

Unruh  and 

Katz 

ch4-co2-h2o 

(52) 

3. 

Reamer  et 

al 

c3h8-c3h6-h2o 

(40) 

The  experimental  results  of  this  work  are  tabulated  in 
Table  I  and  Table  II  and  are  depicted  in  Figure  6.  The 
compositions  indicated  in  this  figure  are  those  of  the 
equilibrium  gas  samples  taken  at  the  point  of  hydrate  formation, 

as  described  in  the  experimental  program. 

The  three  phase  equilibrium  for  a  ternary  system  generates 
a  surface.  This  surface  is  made  of  constant  composition  lines, 
each  representing  HL^G  equilibrium  for  the  particular  composition 
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of  the  gas  phase „  Since  it  is  not  possible  to  manipulate  the 
composition  of  the  gas  phase  at  the  initial  hydrate  formation 
conditions ,  the  points  shown  in  Figure  6  are  scattered.  The 
lines  drawn  in  this  region  are  approximately  located*  All 
these  lines  are  terminated  at  the  quadruple  locus*  It  is 
interesting  to  note  the  occurrence  of  a  maximum  temperature 
of  56*7°F  in  the  quadruple  locus*  There  exist  no  Structure  II 
hydrates  beyond  this  temperature* 

E .  Discussion  of  the  Experimental  Results 

The  general  trend  of  the  results  obtained  in  this  work 
is  in  good  agreement  with  similar  investigations  of  other  research 
workers*  The  addition  of  the  heavy  component  (C^Hg)  lowers  the 
initial  hydrate  formation  pressure  and  temperature  of  the 
light  component  (C02) •  At  the  freezing  point  of  water,  the 
addition  of  about  10%  C.H  changes  hydrate  formation  pressure 

j  O 

of  CC>2  from  175  psia  to  85  psia*  Similar  effects  were 

(9)  (32) 

observed  by  Deaton  and  Frost  and  Otto  during  their 

studies  on  propane-methane-water  and  propylene-methane-water 
systems  respectively* 

Other  ternary  systems  having  phase  behaviour  characteristics 

similar  to  the  C02-H20-CgHg  system  are  H^S-CHCl^-H^O  and 

H0S-CC1  -H  0(56)*  The  quadruple  locii  of  these  systems  look 
1  4  2 

remarkably  alike*  The  points  on  the  four  phase  UL^  G  locus 
are  the  points  of  intersections  of  the  three  phase  branches 
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HL  G,  HL  G ,  HL  L  ,  L  LG  for  a  particular  fixed  gas  composition. 
1  2  12  12 

The  branches  HL^G  and  L^L^G  are  located  below  the  vapor 


pressure  curve  for  the  gas  mixture  of  the  same  composition „ 

f  37  3Q ) 

From  the  available  data'  '  '  on  the  system  CO„-C  H  ,  it 

was  possible  to  draw  dew  point  pressure-temperature  curves  for 
various  compositions.  These  curves  are  shown  in  Figure  7 
along  with  the  quadruple  locus  and  three  phase  HL^G  lines  for 
the  ternary  CC^-C^Hg-H^O  system.  The  three  phase  lines  in 
Figure  7  are  located  approximately  due  to  the  scatter  in  the 
points,  as  shown  in  Figure  6.  By  this  method,  the  shape  of 
the  quadruple  locus  can  be  defined  approximately. 

The  slope  of  a  univariant  phase  equilibrium  branch, 
plotted  on  a  P-T  plane,  is  given  by  the  Glapeyron  equation, 


dP  AH 


dT  T  *  AV 


(15) 


This  thermodynamic  identity  is  valid  irrespective  of  how  the 
compositions  of  the  phases  change  along  the  equilibrium  line. 
The  quantities  AH  and  AV  denote  the  changes  in  the  enthalpy 
and  volume  accompanying  the  phase  transformation  along  the 
univariant  equilibrium  considered  at  constant  P  and  T. 

Approximately,  AV  =  zRT/P  (16) 

Substituting  equation  (16)  in  (15)  and  rearranging, 

d  In  P  _  A  H  (17) 

dT  z  R  T2 

Thus  when  In  P  is  plotted  against  T  on  a  semi-log  graph, 
the  slope  of  this  line  is  given  by  the  right  hand  side  of 
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FIG-7  FOUR  PHASE  -  HDL|L2G- LOCUS  FOR  COj-CjHg-HjO  SYSTEM 
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equation  (17),  namely  by  the  quantity  (AH/zRT  ). 

The  slope  of  the  HL-.G  curve  for  the  C-H0-H  0  system  is 

greater  than  that  of  the  HL^G  curve  for  the  CC^-H-O  system. 

This  is  evident  from  the  fact  that  AH  value  for  carbon  dioxide 

(56) 

hydrates  is  14.4  k-cal/mole  ,  while  the  same  for  propane 
hydrates  is  32.0  k-cal/mole .  This  can  also  be  deduced 
from  the  knowledge  that  the  Afji  value  for  Structure  II  hydrates, 
is  greater  than  that  for  Structure  I  hydrates. 

The  transformation  along  the  four  phase  HL^I^G 
equilibrium  may  be  represented  by, 

LX  +  L2  +  G  - J-  HZI 

The  slope  of  this  univariant  quadruple  phase  equilibrium 
changes  its  direction  at  the  maximum  temperature  point.  A 
tangent  moving  along  this  locus  will  have  maximum  slope  at 
this  maximum  temperature  point  beyond  which  H  phase  does 


not  exist. 
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THEORETICAL  RESULTS 


A.  Potential  Models  for  Non-Spherical  Molecules 

Various  potential  functions  have  been  proposed  to  account 

for  the  non-spherical  shape  of  molecules  like  CC^.  The  parameters 

of  such  a  potential  function  could  not  be  computed  from  first 

principles.  They  are  customarily  obtained  from  some  macroscopic 

property.  Of  these,  the  second  virial  coefficient  itself  has 

been  used  most  conveniently.  Consequently  the  experimental 

information  on  the  second  virial  coefficient  of  various 

molecules  was  collected  to  evaluate  parameters. 

A  considerable  improvement  over  the  Lennard-Jones 

(38) 

potential  has  been  suggested  by  Kihara  who  replaced  the 

point  center  model  for  a  molecule  by  an  impenetrable  core 

whose  dimensions  are  suggested  by  the  geometry  and  inter- 

(7) 

molecular  distances  of  the  molecule.  Corner  proposed  a 

four-center  model  for  long  molecules  like  CO^  and  C^H^. 

When  these  models  were  used  to  predict  the  dissociation  pressures 

of  gas  hydrates,  the  discrepancy  was  quite  pronounced.  It  was, 

(29) 

therefore,  concluded  by  McKoy  that  the  dissociation 

pressures  are  more  sensitive  to  the  parameters  of  the  potential 
functions  than  are  second  virial  coefficient  data.  The 
parameters  for  CO^  and  CgH^  are  tabulated  in  Table  V, 

There  is  a  unique  set  of  parameters  for  a  given 
molecule.  They  could,  therefore,  be  computed  by  fitting 
the  observed  value  of  hydrate  formation  pressure  and  temperature 
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at  a  point  on  the  HL^G  equilibrium  curve.  These  values  of 
(  6 A)  and  (<S)  for  CC^  and  C^Hg  were  used  to  calculate  the 
rest  of  the  pressures  and  temperatures  along  the  three  phase 
(HL^G)  equilibrium.  The  method  of  calculating  these 
potential  parameters  is  described  in  Appendix  D.  The  values 
obtained  by  this  procedure  are. 

Molecule  (  6  A)  /  °K  (6)  ,  A 

C02  148  3.96 

t 

C3H8  278  4,80 

B.  Langmuir  Constants 

The  values  of  Langmuir  constants  by  equation  (3)  are 
dependent  on  the  parameters  t5  and  £A  which  in  turn  are 
characterized  by  the  potential  model  used.  The  parameters  in 
Table  V  were  used  to  calculate  C-^  values  of  C02  and  C^Hg 
molecules.  The  procedure  for  calculation  of  Langmuir  constants 
is  shown  in  Appendix  C„ 

Langmuir  constants  for  C  H  are  plotted  as  a  function  of 

3  8 

temperature  in  Figure  8  and  Figure  9.  It  is  interesting  to 

o 

notice  the  influence  of  change  of  the  value  of  o  from  4.90A 

o 

(Curve  B  of  Figure  9)  to  5.50A  (Curve  B  of  Figure  8),  keeping 
(  BA)  constant  at  228°K.  A  similar  effect  is  brought  about 
due  to  the  change  in  (  6A)  values  at  constant  ( <5  )  value, 
as  shown  by  Curve  A  of  Figure  8  and  Figure  9 .  When  (tf )  is 
increased  by  89%  at  constant  (  6  A)  =  228 °K  the  Langmuir 
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FIG-8  LANGMUIR  CONSTANT  FOR  PROPANE  MOLECULE  IN 
STRUCTURE  II  HYDRATES,  LARGE  CAVITY 
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FIG-9  LANGMUIR  CONSTANT  FOR  PROPANE  MOLECULE  IN 
STRUCTURE  E  HYDRATES,  LARGE  CAVITY 
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constant  increased  by  an  average  of  7.85%  over  the  temperature 
range  of  0-17 °C.  But  when  ( 8 /k)  values  were  raised  by  54% 

o 

at  constant  <5=  4.80  A,  the  C^.  values  increased  by  an  average 
of  0.62%  over  the  same  temperature  range.  Thus  Langmuir 
constants  are  more  sensitive  to  the  ( 6  )  values  than  to  the 
(  6  A)  values  . 

Langmuir  constants  of  the  CO2  molecule  in  Structure  I  are 
shown  in  Figure  10  and  Figure  11  for  various  potential  models. 
For  the  models  reported  in  the  literature,  C^.  values  in  large 
cavities  are  smaller  than  those  in  small  cavities  at  all 
temperatures.  For  the  set  of  values  obtained  in  this  work, 
values  in  large  cavities  are  greater  than  those  in  small 
cavities.  This  can  be  explained  by  equation  (2)  which  relates 
the  probability  of  finding  a  molecule  in  large  or  small 
cavities  to  its  Langmuir  constants.  To  have  yk2  Ykl'  f°r 
a  gas  molecule  in  Structure  I  hydrates,  one  must  have,  . 

The  probability  of  finding  a  molecule  in  the  large  cavities  is 
obviously  greater  than  that  of  finding  the  same  molecule  in  the 
small  cavities  of  Structure  I  hydrates.  The  same  trend  was  also 
observed  by  Marshall^ for  CH^,  A  and  N2  in  Structure  I. 

Langmuir  constants  of  C02  and  C^Hg  are  tabulated  in 
Table  VI  and  Table  VII  respectively. 
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FIG-IO  LANGMUIR  CONSTANTS  FOR  CARBON  DIOXIDE 
MOLECULE  IN  STRUCTURE  I  HYDRATES 
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FIG  II  LANGMUIR  CONSTANT  FOR  CARBON  DIOXIDE  MOLECULE 

IN  STRUCTURE  -  I  HYDRATES 
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C .  Prediction  of  Initial  Hydrate  Formation  Conditions 

Equations  (12)  and  (14)  were  used  to  calculate  initial 

hydrate  formation  conditions,  as  described  in  Appendix  E. 

In  addition  to  the  Langmuir  constants,  values  of  Ap,’  are 

required  for  these  calculations. 

The  chemical  potential  difference  Ap,'  for  HL-^G  equilibrium 

for  Structure  I  and  Structure  II  hydrates  were  calculated  from 

equation  (14)  using  methane  and  sulfur  hexafluoride 

respectively  as  reference  components.  Data  fox  the  methane 

hydrate  were  obtained  from  Shozaburo  et  al^^  and  those  for  SF^ 

were  computed  using  equation  (12)  and  experimental  data  of 
( 48) 

Sortland'  .  These  values  are  tabulated  in  Table  VIII  and 
Table  IX. 

The  initial  hydrate  formation  conditions  predicted  by 
this  method  are  tabulated  in  Table  X  and  are  compared  with 
the  experimental  data  in  Figure  12 . 

D .  Discussion  of  Theoretical  Results 

The  solid  solution  theory  was  limited  by  the  assumptions 
made  during  its  development.  These  assumptions  restricted  the 
size  and  shape  of  the  molecules.  Non-spherical  molecules 
like  C02  will  not  be  free  in  the  entire  cavity.  When  such  a 
molecule  comes  close  to  the  wall  of  its  cage  it  will  have  to 
orient  itself  parallel  to  this  wall.  Furthermore,  some  of  the 
cavities  are  somewhat  oblate  in  shape,  and  thus  the  rotations 
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of  relatively  large  and/or  non-spherical  molecules  may  be 
seriously  hindered.  These  effects  will  predict  higher 
dissociation  pressures  than  those  calculated  by  assuming  free 
rotations . 

(56,  57) 

Platteeuw  and  van  der  Waals  calculated  Langmuir 

constants  of  propane  at  -3°C  using  equation  (2),  in  which 
they  obtained  y  values  for  SFg  molecules  and  dissociation 
pressure  from  the  experimental  data  for  C  H  hydrates.  This 
implies  that  the  equilibrium  composition  of  the  hydrate  is 
independent  of  the  nature  of  the  guest  molecules.  Provided 
that  the  guest  molecules  do  not  distort  the  host  lattice  and 
that  the  pressure  is  low  enough  for  the  variation  of  A|i 
with  pressure  to  be  negligible,  the  Aji  value  will  be  constant 
for  any  molecule,  at  that  temperature  and  for  the  particular 
structure . 

(26) 

Marshall  computed  parameters  ( (f ,  £/k)  for  the  methane 

molecule  by  fitting  different  sets  of  (  & ,  £/k)  values  to  the 
experimental  data  at  more  than  one  point.  The  combination 
(6*/  £A)  which  gave  a  minimum  mean  square  deviation  from  the 
experimental  data  was  chosen.  The  parameters  for  propane 
calculated  by  the  method  proposed  in  this  work  are  in  good 
agreement  with  those  calculated  by  Marshall.  The  Corner 
potential  was  the  one  which  gave  C  values  close  to  the  values 
used  in  this  work. 

To  account  for  the  non-ideal  behaviour  of  CC>2  and  C3HQ, 
fugacities  were  used  in  all  the  calculations.  The  initial 
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hydrate  formation  conditions  are  sensitive  to  the  parameters 

{6)  and  (  6/k)  because  these  parameters  determine  the  Langmuir 

constants  which  at  a  given  temperature  influence  dissociation 

pressures  according  to  equation  (12).  It  was  observed  in  general 

that  the  increase  in  C  values  shifted  the  HL, G  curve  down  on 

k  1 

the  pressure-temperature  plane,  and  vice  versa. 

The  use  of  methane  and  sulfur  hexafluoride  as  reference 
components  proved  very  helpful.  This  approach  assumes  that 
the  guest  molecules  do  not  distort  the  host  lattice.  This 
assumption  was  shown  to  be  a  very  good  approximation  by 
McKoy^29^ . 

The  agreement  between  the  predicted  hydrate  formation 

conditions  and  the  experimental  observations,  as  shown  in 

Figure  12,  is  within  6%.  The  slope  of  the  theoretical  HL^G 

equilibrium  curve  increases  with  pressure  and  temperature. 

This  can  be  explained  by  equation  (17).  It  was  found  that 

compressibility  factor  z  decreases  along  the  HL^G  equilibrium 

curve,  with  the  increase  in  pressure  and  t emperature .  It 

can  be  intuitively  said  that  the  difference  in  the  specific 

heats  of  the  products  and  reactants  is  positive  because  the 

products  (hydrates)  are  solids  while  the  reactants  are  liquid 

and  gaseous.  Thus,  by  Kirchoff's  law,  AH  will  increase  with 

the  increase  in  temperature,  assuming  it  to  be  independent  of 

2 

the  change  in  pressure.  The  change  in  the  T  term, where  T  is 
absolute  temperature,  is  quite  small.  For  these  reasons,  (AH/z) 
will  increase  with  temperature  and  hence  the  slope  of  HL-j^G 
equilibrium  increases  with  the  rise  in  pressure  and  temperature. 


; 

. 

■.  . v  *  >.v 


'.V  bn>  s&a& .r.q  •  ‘ 3 i/ijB : w  ri  7  r.,-  f 


l  rx-  .  rf; 


.  Hi  nre  cijr.p  ai  , -.X&3&H*  c  )  if  loacfr 


CONCLUSIONS 


The  three  phase  (HL-^G)  equilibrium  for  the  binary 
CO2--H2O  system  was  determined  in  the  temperature  range 
32 °F  to  50.3°F.  The  quadruple  point  for  this  system 
was  observed  at  T  =  50.3°F  and  P  =  648  psia. 

The  binary  system  C3Hg-H20  was  studied  in  its  HL-^G 
equilibrium  from  32 °F  to  42.3°F  and  its  four  phase 
(HL2L2G)  point  was  found  at  T  =  42.3°F  and  P  =  80  psia. 

The  three  phase  equilibrium  conditions  for  the  ternary 
C02-C3H8-H20  system  were  determined  at  different 
compositions .  The  quadruple  locus  for  this  ternary 
system  was  investigated.  The  maximum  temperature 
point  was  observed  in  this  quadruple  (HL^L^G)  locus  at 
T  =  56 .7  °F . 

| 

The  parameters  for  the  potential  function  of  CO2  and 
C3H8  computed  in  this  work  were  found  to  be, 

o 

Carbon  dioxide:  (o)  =  3.96  A 

(  6A)  =  148.0  °K 

O 

Propane:  (tf)  =  4.80  A 

(  8 A)  =  278.0  °K 

The  predicted  initial  hydrate  formation  conditions  for 
C02-H20  and  C3H8”H20  systems  in  HL^G  equilibrium  agreed 
satisfactorily  with  the  experimental  data.  Thus,  the  solid 
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solution  theory  was  found  to  be  valid  for  the  three 
phase  (HL^G)  equilibria  of  systems  containing  non- 
spherical  molecules  like  CO  and  C_H0. 

2.  3  o 
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AREAS  FOR  FUTURE  WORK 


To  verify  the  shape  of  the  P-T  diagram  of  the  C02-C3Hg-H20 
system,  some  pressure-composition  envelopes  may  be 
obtained  experimentally.  Such  a  study  may  also  prove  the 
coexistence  of  the  two  solid  phases  H ^  and  H-j-j.  Such 
pressure-composition  data  may  also  be  used  to  test  the 
validity  of  the  solid  solution  theory  for  the  HL^G 
equilibrium. 

To  avoid  the  trial  and  error  procedure  in  calculating 
parameters  a  and  e/k,  attempts  may  be  made  to  fit  the 
data  at  two  points  and  solve  the  simultaneous  equations 
for  two  unknowns. 

With  these  parameters,  Langmuir  constants  may  be  computed 
for  each  hydrate  former.  These  calculations  may  be 
compared  with  the  published  data  on  various  hydrates. 

One  may  theoretically  calculate  the  vapor-solid 
equilibrium  constants  and  pressure-composition 
envelopes  at  desired  temperatures  and  the  pressure- 
temperature  diagrams  for  multi-component  systems. 

With  these  procedures,  one  may  evaluate  the  selectivity 
of  the  hydrate  phase  for  a  particular  component  compared 
to  its  concentration  in  the  equilibrium  gas  phase.  These 
calculations  may  decide  upon  the  possibility  of  using  the 
hydrate  formation  process  for  separation  of  multi-component 


mixtures . 
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APPENDIX  A 


CALIBRATIONS 

A.  Chromatograph  Calibrations 

v 

A  gas  chromatograph,  model  K-2  supplied  by  the  Burrell 

? 

Corporation,  was  used  to  analyze  the  mixtures  of  carbon  dioxide 
and  propane.  Separation  of  both  components  was  satisfactory, 
using  the  following  column  and  conditions? 


Columns  250  cms  long,  packed  with  30-60  #  mesh  Fire  Brick, 
coated  with  Apiezon  M  Grease,  using  Benzene  as  a 
solvent.  The  proportion  used  was  30  gms  of  Fire 
Brick  to  7  gms  of  Apiezon  M  Grease. 

Column  Temperatures  40 °C 


Detector  Currents 


Helium  Rates 
Sample  volumes 


180  mA 


25  ml/min 


Residence  time  for  CO^ s 


0.875  ml 
2  mins 

Residence  time  for  C3HQ  s  3  mins 

Samples  of  known  composition  were  prepared  in  the  gas 
burette  and  analyzed  at  the  above  conditions .  Peak  height  of 
each  component  was  plotted  against  its  composition.  Unknown 
samples  were  then  tested  and  their  compositions  normalized 
on  a  prorata  basis.  The  calibration  curve  is  shown  in  Figure 
13.  This  curve  is  on  the  basis  of  attenuation  of  50.  Accuracy 

_L, 

of  this  calibration  is  within  -  3%.  The  analysis  of  lean 
samples  is  not  very  accurate  because  the  peak  height  of  the  lean 
component  was  not  very  reliable. 
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FIG-13  CHROMATOGRAPH  CALIBRATION  FOR  C02~  C3H8  MIXTURES 
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B .  Thermocouple  calibration 

The  thermocouple  was  standardized  at  the  ice  point  and 
room  temperature.  The  calibration  curve  after  corrections  is 
shown  in  Figure  14. 
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APPENDIX  B 


Structural  Details  of  Hydrate  Crystals 

Structure  I  -  The  unit  cell  of  Structure  I  contains  46 
molecules  of  water  which  form  two  types  of  cavities.  The  small 
cavities  are  located  at  the  vertices  and  center  of  the  unit 
cell.  The  small  cavity  is  formed  by  a  pentagonal  dodecahedral 
array  of  20  water  molecules.  The  remaining  six  water  molecules 

I 

in  the  unit  cell  form  bridges  between  these  dodecahedra  in 
such  a  way  that  a  second  type  of  cavity,  a  large  cavity,  is 
formed.  The  large  cavity  is  a  tetrakaidecahedron  having  two 
opposite  hexagonal  faces  and  twelve  pentagonal  faces.  Thus, 
Structure  I  has  the  following  constants; 

Water  molecules  =  =  46 

Small  cavities  =  n^  =  2 

Large  cavities  =  ^  =  6 

Number  of  oxygen  atoms  surrounding  small  cavity  =  =  20 

Number  of  oxygen  atoms  surrounding  large  cavity  =  z2  =  24 

Average  distance  of  oxygen  atoms  f:gom  the  center  of  small 
cavity  =  a^  =  3.9  5  A 


Average  distance  of  oxygen  atoms  f£om  the  center  of 
cavity  =  a2  =  4.30  A 

Number  of  small  cavities  per  mole  of  water  =  *\) 


Number  of  large  cavities  per  mole  of  water  =  \J 


Free  diameter  of  small  cavities  —  5.1  A 

o 

Free  diameter  of  large  cavities  =  5.8  A 


large 
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Chemical  potential  difference 
at  0°C 

Volume  difference  =  AV  = 

=  AV’  = 

Structure  II  -  The  smaller  cavities  in  this  structure  are 
distorted  pentagon  dodecahedra  and  the  larger  cavities  are 
almost  spherical.  Various  constants  for  this  structure  ares 
Water  molecules  =  Nq  =  136 

Small  cavities  =  n^  =  16 

Large  cavities  =  n2  =  8 

z1  =  20,  =  3.91  A,  =  2/17 

7*2  -  28,  a 2  =  4.73  A,  =  1/17 

O 

Free  diameter  of  small  cavity  =  5.0  A 


Free  diameter  of  large  cavity  =  6.7  A 


A(j,  = 

a 

-  = 

190  cal/mole  at  0°C 

AV  = 

vP 

-  va  = 

3.4  ml /mole 

§ 

11 

vP 

il 

a 

> 

1 

5.0  ml/mole 

The  above  values  are  quoted  from  reference  56. 


=  An,  =  hP  -  na  =  167  cal/mole 

Q  Q 

-  Va  =  3.0  ml/mole 

-  VL  =  4.6  ml/mole 
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APPENDIX  C 


Calculations  of  Langmuir  Constants: 


Langmuir's  constant  is  defined  by  the  equation, 


C.  • 
ki 


where, 


2tt  a?  g  . 
1  ki 


k  T 


exp 


wki(°) 


kT 


1/4 


g,  . 

ki 


where , 


exp 


zi6 


k  T 


(-- 


(y)  2m( y) 

+  - : 


4 


a 


•  y 


1/2 


dy 


1  ( y )  =  j^(y4  +  12y 3  +  25. 2y2  +  12y  +1)  /  (1  - 


y) 


-  l 


m 


(y)  =  Q  (i  +  y)  /  (i  -  y)4^] 


-  l 


ai  =  {a±  /  6  )3  /V^2 

6  -  (  tfQ  +  /  2 


6A 


wkj(0) 


k  T 


„  1/2  1/2 

(  £A)q  (  6A)k 


z±  6 


k  T 


4  2 

ai  ai 


To  calculate  the  g  ,  "free  volume",  numerical  integration 

ki 

by  using  Simpson's  Rule  was  carried  out  on  a  digital  computer 
IBM  7040.  According  to  Simpson's  Rule, 


T  X 


\ 


r*<,  - 1)  \  (,  *  i)  i 


£  \  (  ,b  .&  ) 


A^ 


T  tf 
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Y0  +  4  (Yl  +  Y3  +  *  *  *  *  +  Yn-l)  + 

2  (y2  +  y4  +  ...  .  yn_2)  +  yn 

The  area  under  the  curve  is  divided  into  n  even  number  of 
slices,  each  of  a  width  equal  to  h.  In  the  present  case, 
a  =  0,  b  =  1/4,  yQ  =  0.  The  integrand  decays  as  y  increases 
and  hence,  it  was  found  that  y^  =  0 .  With  these  simplifications, 
Simpson's  Rule  becomes, 


b 


f(y)  *  d  y  =  — 

3 


f  (y)  *  d  y  = 


2h 

3 


'  2  (Yj_  +  Y3  +  •  • 


•  •  +  Y  n  )  + 

Jn-1 


(y2  +  y4  +  *  ’  ’  •  +  ^ 

In  this  work  the  area  was  divided  into  1024  parts  when  the 
values  of  g^.  converged  in  the  fifth  decimal  place. 

The  Langmuir  constant  is  made  up  of  three  terms,  namely 

(i) 

(ii)  g.  .  and 

(iii)  exp 


(2ir  /  kT) 

gki  r-  ,  ^ 

 Wki(0) 


k  T 


Term  (i)  is  a  volume  term  and  is  a  constant  for  a  given 
structure  and  cavity  type.  Term  (ii)  accounts  for  the  free 
volume  in  a  cavity.  Term  (iii)  is  an  exponential  function  and 
hence  is  sensitive  to  the  values  of  parameters  6  and  S/k.  As 
the  value  of  6  increases,  the  Langmuir  constant  increases.  A 
similar  effect  is  observed  due  to  an  increase  in  the  ( 6  /k) 
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values.  These  are  evident  in  Figures  8  to  11.  In  all  cases 
Langmuir  constants  decrease  with  increasing  temperature. 
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APPENDIX  D 


Calculation  of  Parameters  <5  and  S/k 

For  the  solid  solution  theory  to  be  exact,  there  should 
exist  a  unique  set  of  values  of  6/k  and  ^of  a  gas  molecule 
for  all  its  data  on  hydrates.  Once  these  values  are  determined, 
Langmuir  constants  could  be  calculated  as  a  function  of  temperature 

The  reported  values  of  these  parameters  show  large 
variations  as  shown  in  Table  V,  In  order  to  fit  these  constants 
to  the  hydrate  data,  an  attempt  was  made  to  force  fit  the 
experimental  and  calculated  values  of  dissociation  pressure  at 
0°C.  As  there  are  two  unknowns,  £/k  and  one  of  the  two 

was  assumed,  and  the  other  one  adjusted  by  the  iteratative 
method  until  the  observed  and  calculated  values  of  the 
dissociation  pressure  agreed.  In  this  work  6  was  assumed  for 
each  gas  molecule  and  (  S’/k )  adjusted  until  the  two  values  of 
dissociation  pressure  matched. 

The  solubility  data  of  carbon  dioxide  (19,  58,  59) 

and  Propane  (21,  24)  were  extrapolated  into  the  hydrate 


region  by  Othmer ' s  Plot 


(33) 


.  Sulfur  hexafluoride  was  found 


(14)  to  be  negligibly  soluble  in  water.  Data  for  methane 

(46  47) 

hydrates  were  obtained  from  Shozaburo  ' 


Carbon  dioxide,  (^)  and  (  6/k) 

Pure  carbon  dioxide  grows  in  Structure  I  and  occupies 

its  large  and  small  cavities. 
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From  equations  (12)  and  (14),  we  have, 


A  iul  • 


A  (i 


RT 


^1  ln  U  + 


A|jl  1  + 

o 


AV  '  (P 


Po’ 


+  \)  In  ( 1  +  C  p  ) 

v  2  k2 


-i-  In  x 


Q 
(12) 
(14) 


Selecting  methane  as  a  reference  component,  values  of 

(46,  47) 

Afj/  and  PQ  were  obtained  from  Kobayashi  et  a.1  .  In 

order  to  account  for  the  nonideality  of  carbon  dioxide, 
fugacities  (p^)  were  used  in  all  the  calculations.  From 
experimental  observations  it  was  found  that  P  =  172.0  psia  at 

o 

0°C.  Assuming  5  =  3.96  A,  Ckl  and  C^2  were  calculated  for 
every  iteratative  value  of  (  6  /k) .  Knowing  the  fugacity 
at  P  =  172.0  psia  and  T  =  0°C,  and  the  solubility  (x^),  values 
of  Ajjl  '  from  equation  (12)  were  computed  on  the  IBM  7040 
computer  for  each  iteration.  This  process  of  iteration  was 
continued  until  A|j, '  values  from  equations  (12)  and  (14) 
matched.  The  parameters  so  determined  for  carbon  dioxide  were 

^  O 

(O)  =  3.96  A,  (  S  A)  =  148°K.  The  Langmuir  constants  using 
these  parameters  are  shown  in  Figure  11. 

Propane,  (  6  )  and  (  6  A) 

Propane  forms  Structure  II  hydrates  and  enters  the  large 


cavities  only. 
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From 

*■ 

equations 

(12)  and  (14),  we  have 

ri 

<1 

=  RT 

V2  in  (1 

+  Cp2pp)  +  ln  XQ 

Ap,' 

=  And, 

+  A  V' 

(P  -  PJ 

(12) 

(14) 


Sulfur  hexafluoride  was  chosen  as  a  reference  component 

for  Structure  II  hydrates.  Pressure-temperature  experimental 

(48) 

data  for  SF^  hydrates,  as  reported  by  Sortland  were  used. 

Langmuir  constants  for  the  SF^  molecule  are  given  in  Table  IX. 
Using  these  values,  Ap,^  was  calculated  from  equation  (12)  when 
Xq  =  1.0  because  SF^  is  negligibly  soluble  in  water.  Assuming 

o 

< 5  =  4.80  A  for  the  propane  molecule,  (  £ /k)  was  determined 

by  the  iteratitive  procedure,  similar  to  the  one  described  for 
carbon  dioxide.  Values  of  these  parameters  so  obtained  were 
6  =  4.80  A  and  (  6  A)  =  278  °K. 
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APPENDIX  E 


Prediction  of  Initial  Hydrate  Formation  Conditions 

With  the  best  fitted  values  of  parameters  €>  ,  6  / k 

for  carbon  dioxide  and  propane,  as  described  in  Appendix  D, 
Langmuir  constants  were  calculated  as  a  function  of  temperature. 
(See  Figures  9  and  11) . 

At  a  selected  temperature,  values  of  Ap,'  can  be  calculated 
for  each  structure,  at  various  total  pressures  P,  from  equation 
(12).  At  the  same  temperature,  values  of  Langmuir  constants 
are  known.  From  solubility  data,  Xq  is  determined  at  this 
temperature  and  various  total  pressures  P.  An  iteratative 
procedure  was  set  up  for  various  values  of  the  fugacity  of 
the  pure  component  at  the  same  temperature  and  various  total 
pressures  P,  to  calculate  Ap,’  from  equation  (12).  The 
iterations  were  continued  until  the  values  of  the  chemical 
potential  difference  calculated  from  equations  (12)  and  (14) 
matched . 

Results  of  such  calculations  are  shown  in  Table  VI, 
for  and  C3H8~H2°  systems. 
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TABLE  I 


Experimental  Data  on  Initial  Hydrate  Formation  for 
C02~C2Hg-H20  System  in  Three  Phase  -  HL-^G  Equilibrium 


Mole  %  CO 2  in 

Pressure, 

Temperature, 

Mole  %  CO 

Feed  Gas 

psia 

op 

G-phase 

100 

200 

33.4 

100 

226 

35.7 

100 

255 

37 .3 

100 

351 

42.4 

100 

454 

45.6 

100 

557 

48.0 

100 

95 

461 

51.0 

92.0 

317 

47  .4 

92.0 

182 

41.3 

90.0 

120 

36.0 

85.0 

98 

33.4 

92.0 

93 

619 

53.0 

94.5 

440 

50.6 

91.0 

2  57 

45.0 

91.0 

211 

42.3 

91.0 

167 

37  .6 

94.0 

118 

33.2 

93.0 

91 

490 

55.4 

84.0 

228 

46.0 

87.0 

175 

43.3 

86.0 

110 

38.0 

85.0 

75 

33.2 

87.0 

82 

278 

51.1 

75.0 

142 

44.7 

75.0 

95 

39.6 

76.0 

52 

33.5 

79.0 

80 

189 

47.6 

74.0 

103 

40.8 

77.0 

in 
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Mole  %  CO2  in 

Feed  Gas 

Table  I  - 

Pressure , 

psia 

( continued) 

Temperature , 

op 

Mole  %  CO 

G-phase 

64 

240 

51.0 

58 

155 

46.3 

52 

100 

41.8 

52 .5 

60 

36.5 

57  .8 

46 

115 

43.8 

40 

109 

43.6 

37 

84 

41.2 

35 

63 

37  .8 

39 

47 

35.0 

40 

29 

86 

42.6 

18 

73 

41.1 

17 

44 

35.7 

28 

15 

93 

42  .7 

17 

69 

40.3 

16 

49 

37 .1 

14 

0 

66 

40.4 

0 

48 

37  .8 

0 

35 

35.0 

0 

30 

34.0 

0 

in 
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TABLE  II 


Experimental  Data  for  Four  Phase  HL^L2G 
Equilibrium  for  CC^-CgHg-F^O  System 


Mole  %  C02  in 

Pressure, 

Temperature , 

Mole  %  < 

Feed  Gas 

 psia 

°F 

G  Ph 

100 

648 

50.3 

100 

95 

632 

55.2 

91 

91 

580 

56.7 

87 

82 

435 

55.5 

76 

80 

462 

56.4 

82 

64 

300 

52  .9 

60 

46 

169 

48.6 

49 

29 

139 

46.6 

37 

0 

80 

42.3 

0 

TABLE  III 


Experimental  Data  of  Other  Workers  for 
C^-HpO  System  in  HL^G  Equilibrium 


Authors 


Pressure , 
psia 


Temperature, 

op 


Unruh  and 


(52) 

Katz 


653 

49  .6 

535 

47  .7 

468 

45.9 

37  5 

42  .8 

296 

39.3 

Maximum  temperature  for  HL^G  equilibrium  =  50.36°F 


(53) 


179.34 

186.69 

245.49 

263.13 

320.46 

342.51 

383.67 

651.21 


32 .00 

32 .86 

36.86 
38.48 
41.54 
42  .98 
44.42 
50.00 


in 


Villard 


ri 


■  .  : 


Itosi  txA 

(£6) 
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TABLE  IV 


Experimental  Data  of  Other  Authors  for 
C^Hg-I^O  System  in  the  Three  Phase  HL^G  Equilibrium 


Authors 

Pressure, 

Temper  at- 

psia 

op 

82.0 

42  .25 

79.5 

42.07 

69.5 

40.91 

Miller  and 

68.5 

40.82 

(34) 

66.5 

40.46 

Strong 

57 .0 

39.20 

56.5 

39.11 

53.0 

38 . 66 

50.0 

38.30 

(99.9%  pure  propane) 

46.0 

37  .49 

42.5 

36.86 

32 .0 

34.52 

27.5 

33.26 

27.0 

33.08 

25.5 

32  .63 

25.0 

32.45 

24.0 

32.09 

58*24 

38.8 

Deaton  and 

42  .22 

36.3 

29.72 

33.3 

Frost  34^ 

27  .97 

33.1 

27  .84 

32 .8 

(99.8%  pure  propane) 

26.28 

32.4 

Reamer,  Selleck  and 

60.0 

39.3 

44.2 

36.6 

c  (40) 

Sage 

34.9 

34.1 

tosUuA  lediO  lo  £.i&Q  I&Ja9m±x&qx& 


vAIPCl^l 


I.  S6 
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TABLE  V 

Values  of  Force  Constant  £/k  and 
Molecular  Parameter  ($ 


o 


Molecule 

(  S  A) ,  °K 

(<J),  A 

Source 

Carbon 

380.2 

3.95 

(38) 

dioxide 

309.0 

3.36 

(29) 

198.0 

3.57 

(18) 

148.0 

3.96 

This  work 

Propane 

515.14 

4.80 

(38) 

228.00 

5.50 

(  7) 

228.00 

4.90 

(18) 

278.00 

4.80 

This  work 

Sulfur 

200.9 

5.51 

(18) 

Hexafluoride 

Methane 

150.55 

3 .66 

(46) 

•  :  v  \j. 


.OcfifiO 


9fl  <qcr:r  % 


. 
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TABLE  VI 


Langmuir  Constants  for  CO^  in  Structure  I 
Hydrates,  Using  (5)  =  3.96  A,  (  8/k)  =  148  °K 


Ck,  small  cavity, 
Temperature,  °C  _ per  atm. _ 


Ck,  large  cavity, 
per  atm. 


0 

0.3934 

0.5095 

1 

0.3796 

0.4933 

2 

0.3665 

0.4778 

3 

0.3538 

0.4629 

4 

0.3417 

0.4486 

5 

0.3301 

0.4347 

6 

0.3190 

0.4214 

7 

0.3083 

0.4086 

8 

0.2980 

0.3963 

9 

0.2882 

0.3844 

10 

0.2787 

0.37  30 

11 

0.2696 

0.3620 

12 

0.2609 

0.3513 

13 

0.2525 

0.3411 

14 

0.2445 

0.3312 

15 

0.2367 

0.3217 

16 

0.2292 

0.3125 

17 

0.2221 

0.3036 
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TABLE  VII 

Langmuir  Constants  for  C_H  in  Structure  II 

3  8 

Hydrates,  Large  Cavity,  Using  (6  )  =  4.80  A,  (  8 /k)  =  278  °K 


Temperature,  °C 

C  ,  per  atmosphere 
k 

0 

226.99 

1 

214.76 

2 

203.26 

3 

192.46 

4 

182.30 

5 

172  .75 

6 

163.76 

7 

155.29 

8 

147 .32 

9 

139.81 

10 

132.73 

11 

126.05 

12 

119.75 

13 

113.81 

14 

108.20 

15 

97  .89 

16 

93.16 

17 

88.69 

,  •  ■ 
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TABLE  VIII 


Chemical  Potential  of  Methane 
Hydrate  in  Structure  I 


p 

Total 

f 

T 

CM1 

c 

^M2 

Air ' 

cal  per 

Press . , 

Fuqacity , 

Temp . , 

small  cavity, 

larqe  cavity, 

mole 

26.0 

24.5 

273.16 

0.2028 

0.2175 

172.3 

68.0 

59  .2 

282  .90 

0.1525 

0.1690 

231.1 

136.0 

106.2 

289.00 

0.1285 

0.1450 

275.1 

204.1 

148 . 3 

292.30 

0.1173 

0.1332 

301.4 

272.1 

191.0 

294.40 

0.1099 

0.1257 

322 .0 

T,  °K 

Chemical 

Hydrate 

P,  psia 

TABLE  IX 

Potential  of  SF_ 

6 

in  Structure  II 

C  ,  large  cavity,  atm  ^ 

A|r',  cal 
per  mole 

273.16 

11.82 

470.97 

189.64 

275.16 

17.19 

415.31 

199.00 

277.16 

26.18 

366.90 

210.05 

279.16 

41.30 

324.71 

222 .43 

281.16 

66.68 

287.87 

235.80 

283.16 

108.96 

255.65 

249 .79 

285.16 

178.15 

227.41 

264.04 

287.16 

288.07 

202.62 

278.14 

289.16 

455.48 

180.82 

291.71 

291.16 

696.17 

161.61 

304.35 

,  .t.B 

r.*o$: 

<  e.;  ' 

i .  sn 

* 

*•  %T 

dl.C  x 

cl, TVS 

dr  £*£ 

ai.m 

s 

sk ae*> 

cii.ies 
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TABLE  X 


Calculations  of  Initial 
Hydrate  Formation  Conditions 


System 

P,  psia 

T,  °F 

C  H  -H  0 

3  8  2 

25.95 

32 .0 

31.21 

33.8 

37  .80 

35.6 

46.76 

37.4 

58.70 

39.2 

74.28 

41.0 

CVH2° 

191.1 

32 .0 

237  .2 

35.6 

293.5 

39.2 

373.4 

42.8 

496.1 

46 . 4 

650.0 

50.0 

ov.ea 
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NOMENCLATURE 

a 

o1 

A 

Cell  radius  of  i  type  cavity 

Angstrom  unit 

A,  B 

Hydrate  forming  components 

C 

Number  of  components 

CL  • 

ki 

Langmuir  constant  of  k  molecule  in  i  type  cavity 

f 

Fugacity 

F 

Number  of  degrees  of  freedom  or  free  energy  per  mole 

Pki 

Free  volume  of  k  molecule  in  i  type  cavity 

G 

Gas  phase 

AH 

Enthalpy  difference  per  mole 

H 

Hydrate  phase  or  enthalpy  per  mole 

h  ,  h  . 

Ji  Ki 

Partition  function  of  J  (or  k)  molecule  in  i  type 
cavity 

I 

Ice  phase 

k 

Boltzmann's  constant  or  conversion  factor 

L 

1 

Water  rich  liquid  phase 

L2 

Water  lean  liquid  phase 

N 

Number  of  molecules 

n 

Number  of  moles,  number  of  cavities 

P 

Pressure  or  number  of  phases 

P 

o 

Total  pressure  of  a  reference  component 

Pj'  Pk 

Partial  pressure  of  J  (or  k)  component,  fugacity 
when  non-ideal  behaviour 

Q 


R 


Lattice  former,  water 

Gas  constant,  hydrate  former 


t.in~  non  c  c 


A 


e  Jnnnoqiroo  5o  edmuft 


D 


e-.  .  .".  '■ 

a-.-  c  o©  :fq  >yl 


HA 


Y^xvso 


iji1  ‘it, 


i  *••  fq  ‘Dl 


Vvfic  -  o  n  ,  -■  ■  '  -  ■  --s:?  \  on 


•  f  7  nr  • 


c  xup  I  nr>  o  X  •  «vr 

•  Tr  :> .  ->rr;  .  o  ...sc  m/H 


q 


•  d  ,  }?:  t  i.cf) 


Radial  position 
Absolute  temperature 
Volume 

Molar  volume  difference  between  (3  and  L  forms 

Mole  fraction  of  gas  k  in  equilibrium  with  H  and 
phases 

Mole  fraction  of  water  in  phase 

Probability  of  finding  k  molecule  in  i  type  cavity 

Concentration  of  hydrate  former  k  in  H  phase 

Number  of  oxygen  atoms  surrounding  a  cavity  or 
compressibility 


Stable  modification  of  lattice  former  (ice  or  water) 
or  dimensionless  parameter  of  the  cell  geometry 

Metastable  modification  of  lattice  former  (empty 
lattice ) 

Absolute  activity  of  molecule  k 
Cavities  of  type  i  per  mole  of  Q 
Molecular  parameter  ,  force  constant 
Chemical  potential  difference  per  mole 
Molecular  partition  function 

SUBSCRIPTS 


Gaseous  state 
Type  of  cavity 
Any  component 

Carbon  dioxide,  any  component 


.>  •  C  :  Ir  >  1 


O  xoi  .i  •  r  rd'  i--s; i 


t  JO  nu?. 


Liquid  state 
Methane 

Reference  component 
Propane,  constant  pressure 
Lattice  former,  water  molecules 
Solid  state 
Constant  temperature 

SUPERSCRIPTS 


Crystalline  stable  form  of  lattice  former 

Metastable  form  of  lattice  former 

Liquid  stable  form  of  lattice  former 

Chemical  potential  assuming  ideal  solution 

Chemical  potential  difference  assuming  ideal  solution 
behaviour 


.  ' 
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